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Abstract. We consider the Euler equations governing relativistic compress- 
ible fluids evolving in the Minkowski spacetime with several spatial variables. 
We propose a new symmetrization which makes sense for solutions containing 
vacuum states and, for instance, applies to the case of compactly supported 
solutions, which are important to model star dynamics. Then, relying on these 
symmetrization and assuming that the velocity does not exceed some thresh- 
old and remains bounded away from the light speed, we deduce a local-in-time 
existence result for solutions containing vacuum states. We also observe that 
the support of compactly supported solutions does not expand as time evolves. 
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1. Introduction and main result 

The dynamics of relativistic compressible fluids evolving in the Minkowski space- 
time with n spatial variables is governed by the Euler equations (for instance [1]) 



(1.1) 



k—l 



Here, p and u = (wj)i<j<„ denote the mass density and the (n-dimensional) ve- 
locity vector of the fluid and functions of the variables (i, x) G M_|_ x R", while the 
parameter 1/e represents the light speed and 5jk denotes the Kronecker symbol. 
The range of physical interest for the unknown {p, u) is defined by 

n 

(1.2) p>0, |u|2:=^u2<e^2^ 



while the pressure p = p{p) is assumed to satisfy 
(1.3) 0<p'[p)<e-\ 

Under these conditions, it can be checked that the system of conservation laws (|1.1[) 
is symmetric hyperbolic as long as vacuum is avoided, i.e. under the restriction 
p > 0. That is, it can be written in the symmetric hyperbolic in the so-called 
entropy variables (Makino and Ukai [HlG]) 



dtW + y Ak{W)d^,W = 0, 



k=l 



where W :— {p,u) G M-|_ x Si/^ (the set being the open ball with radius 1/e) 
and such that, for every unit vector v = [vk) G M", the matrix X]fe=i ^kAkiW) 
admits real eigenvalues and a basis of eigenvectors. It is also established in 
the initial value problem with non-vacuum initial data, i.e. 

(1.4) {p,u){Q,-)^{-p,u), 

when the initial data (p, w) G M+ x Bi/^ is bounded away from vacuum, admits a 
local-in-time solution. 

In the present paper, we are interested in the symmetrization and the local-in- 
time existence for the relativistic fluid equations (jl.ip when the initial data (p, u) 
take arbitrary values in M+ x Bi/^ and are allowed to contain vacuum states. Some 
partial but pioneering results were obtained on this problem by Rendall [5] and 
Guo and Tahvildar-Zadeh [2] . For an overview of the standard theory, we refer to 
the relevant the chapter on the Euler equations in Choquet-Bruhat's book [1]. 

We emphasize that compactly supported solutions are important in the applica- 
tions, for instance to model the dynamics of stars. However, the current existence 
theory does not cover this situation; indeed, the transformation proposed in [51 [5] 
does not apply for our purpose since the coefficients Ak{W) therein blow-up near 
the vacuum. On the other hand, when e = 0, the system (|l.ll) reduces to the non- 
relativistic Euler equations, for which local existence of solutions, even solutions 
containing vacuum states, was established by Makino, Ukai, and Kawashima [7]. 
The objective of the present paper is precisely to provide a suitable generalization 
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of the theory in ^ to encompass relativistic fluids. This will be achieved by in- 
troducing yet another symmetrization which significantly differs from previously 
proposed ones. 

An outline of this paper is as follows. In Section [2l we begin our investigation of 
the Euler equations and derive a first version of our symmetrization, which is based 
on using as main unknowns "generalized Riemann invariants" and a "normalized 
velocity" . Then, in Section [3l we recall basic material on Lorentz transformations 
and establish a technical lemma. Next, in Section 3] we are in a position to establish 
the existence result of this paper. The main technical difficulty is to check a positive- 
definiteness property for the symmetric system. Finally, in Section [5] we conclude 
with some remark about the support of solutions. 



2. Symmetrization of the relativistic Euler equations 

Modified mass and velocity variables. We are going to define new variables 
defined by nonlinear transformations of the mass density and the norm of the 
velocity vector, which allow us to put the relativistic Euler equations in a symmetric 
form. To begin with, we introduce the modified mass density variable w by 

w = w{p} / -r^rfs, 

(2.1) Jo 

c(p) Vp'ip)^ lip) ■= P + <^'^Pip), 

where c{p) :— \/p'{p) represents the sound speed in the fluid. 

From now on, we assume that w{p) deflned above is finite. This is the case if, 
near the vacuum, the equation of state is asymptotic to the one of polytropic perfect 
fluids, i.e. p{p) k with 7 > 1 and fc > 0. In the special case p{p) ~ k p^ and 
when e is taken to tend to 0, then the function w{p) approaches p(T~i)/2 (up to a 
multiplicative constant), which precisely coincides with the function introduced in 
[7j in the non-relativistic case. 

In addition, based on the norm \u\ of the velocity vector u, we define the modified 
velocity scalar 

„x /I IX 1 , / 1 + e 

(2.2) V = v{\u\) := — In 



2e Vl-e|M| 

We also introduce the n-dimensional, normalized velocity vector and the associated 
projection operator 

u 

(2.3) , E(u) := I - u(E)u, 

\u\ 

respectively, where / denotes the n x n identity matrix. Observe that E{u) is 
singular as a function of u, when u is close to origin; however, the map |up E{u) is 
actually smooth. 

In the rest of this section, (p, u) denotes a given smooth solution to (|l.ip . 



Proposition 2.1 (Formulation in terms of the modified mass and velocity vari- 
ables) . The relativistic Euler equations are equivalent to the following system in the 
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variables {w, v, u): 

(1 - t^\u\^c{pf) dtw + (1 - c(p)2) u ■ 

+ c{p) (1 - e^lup) u-Vv + c{p) |u|V • u = 0, 
(1 - e^\u\''cipf) dtv + c{p) (1 - e^lup) u ■ Vw 

+ (1 - c{pf) u-Vv~ e^c{pf |upV -5 = 0, 
(1 - eVl^)"^ W\ {dtu + M • Vu) + c{p) E{u) Vw = 0. 
Proof. Step 1. The first equation in p.ip takes the form 



2 



dtp + q{p)-7-^ o I u-dtu-e p [p) dtp 



l_e2|y|2 ■^V^'(l_g2|y|2)2 
+ l„g2|^|2 + g(^) (l_g2|^|2)2 " ' + l_g2[^|2 ^-^'"^ = 

or equivalently 

\ + eM\ci'{p)-\) ^ , ^ 

T^T^H^ ^*^ + T^T>F 

By muhiplying this equation by w' {p) we find 

To rewrite the above equation in a more convenient form, we observe that 

dv 1 



d\u\ l-e2|„|2' 

and so, after further multipHcation by (1 — the equation for the modified 

mass density reads 

(1 + e^\u\^c{pf\ dtW + (1 + c{pf) u ■ Vw 

(2 4) 

+ 2t^c{p) \u\ {dtv + u-Vv) + c{p) V • u = 0. 



Step 2. Next, we expand the second equation in (|f .ip and obtain 

+ ^V-( '^^^l ^ ^)+ '^^^i u-Vu+p'{p)Vp = 0, 

which, after using the mass equation, yields 

Y^^j^{dtu + u-S/u)+p'{p) (e^dtpu + Vp) =0. 

Muhiplying by l/q{p) we arrive at an equation for the velocity vector 
(2.5) (1 - \u\^y^ {dtU + u • Vu) + c{p) {e^dtw u + Vw) = 0. 
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We now multiply ()2.5|) by the vector u itself, and obtain 

(1 -e^ lup)-! (at|up + M- V|wp) + 2c{p) {e^dtw \u\^ +u-Vw) = 0, 
which, after a further multiplication by (2|m|)^^, becomes 
(2.6) dtV + u-Vv + c{p) {e'^\u\dtw + u ■ Vw) = 0. 

Step 3. To derive the equation for u we multiply (|2.5p by the projection matrix 
E{u) and obtain 

(1 - e^lup)"^ E{u){dtu + u ■ Vm) + c{p)E{u){e^ dtwu + \7w) ^ 0. 
In view of the identities 

E{u)u = 0, 

E{u) dfU — dtu — udt\u\ — \u\dtu, 

E{u){u ■ Vw) = \u\u ■ Vu, 

we arrive at the (third) equation for the normalized velocity u, as stated in the 
proposition. 

Finally, we are in a position to return to the equation (|2.4p and, by plugging 
(EH) in dMl, find 

(1 - e^\u\^c{pf) dtw + (1 - c{pY) u-\Iw + c{p) V • u = 0. 

Hence, observing that 

V • u = (1 - e^l^n u • Vw + |u|V • u, 

we obtain the desired equation for w, as stated in the proposition. In turn, we can 
also put (|2.5p in the form stated in the proposition for the function v. □ 

A symmetric hyperbolic formulation. At this juncture, it may be interesting 
to consider the one-dimensional case n — I. Considering the result in Proposi- 
tion l2.1l and setting '■— Vu, etc. and then observing that (in the one-dimensional 
case) E = and u = 1, we obtain the following form of the relativistic Euler equa- 
tions in the variables (w, v) 

(1 - e\^cip)^) dtW + (1 - c(p)2) uu), -f c(p) (1 - eW) - 0, 

(1 - e^U^cipf) dtV + C{p) (1 - £^^2) + (1 _ £2 c{pf)uv., = 0, 

which, obviously, is a symmetric hyperbolic system. In particular, it is obvious that 
the coefficient 1 — e^v?'c{p)'^ remains bounded away from zero, provided the sound 
speed or the fluid velocity scalar (or both) remain bounded away from the light 
speed. 

To derive our symmetrization in general dimension, we need two additional ob- 
servations. 

• First, since u has unit norm we can write 

V • u = V • w - u • = ix{E{u) Vu) , 

where "tr" denotes the trace of a matrix. This allows us to rewrite the last 
term of the w-equation in Proposition [^711 in the form 

c(p) lu|V • u = c(p) |u|tr(£;(u) Vw). 
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Interestingly enough, this term can now be viewed as the "symmetric coun- 
terpart" of the term 

|u| c{p) Eiu) Vw 
aheady contained in the M-equation in Proposition [^TT] 
• At this stage, only one term poses some problem if we are to reach the 
desired symmetric form, that is, the term 

in the u-equation of Proposition 12.11 To compensate for this term, one 
would need to have the term — e^c(p)^ |up Vu in the w-equation, but it does 
not appear that a direct transformation could achieve this. So, we introduce 
a further transformation based based on still some new unknowns: 

(2.7) z± := V ± w, 

which we will refer to as the generalized Riemann invariant variables. Ac- 
cording to Proposition [2T1I we have the equations 

(1 - e^\ufc{pf) dtz± + (1 - c{pf){\u\ ± c{p)) u ■ Vz± 

± (1 =F e^c(/9)|u|) c{p) |u| V ■ u = 0, 

(1 -e^lup)-! \u\'^{dtu + u-S/u) + '^\u\c{p)E{u) (Vz+ - Vz_) = 0. 

Consequently, by combining together the above observations we arrive at the 
main conclusion of this section: 

Proposition 2.2 (Symmetric form of the Euler equations). In terms oj the gener- 
alized Riemann invariant variables (z_|-,z_) and the normalized velocity u defined 
in ()2.ip . (|2.2p . (|2.7p . the relativistic Euler equations take the following symmetric 
form 

(2.8) ' ^ ^■^^'^^>'"""'+ ' l^e\{p)\u 



(1 + e2|u|c(p)) dtz+ + f — + c(p)) u ■ Vz+ 

+ c{p)\u\tr{E{u) Vu) = 0, 



,^ 0, (1 - ^Mc{p)) dtz^ + , ^ , /(^ I (1^1 ~ c[p)) u ■ Vz_ 

(2.9) l + e''c{p)\u\ 

- c{p) \u\ tr{E{u) \/u) = 0, 

(2.10) ^J-p. idtu + u-Vu)+ c{p) \u\ E{u) Vz+ - c{p) \u\ E{u) Vz. = 0, 

1 — e^\uY 

where the unknowns z± are real-valued and u is a unit vector, \u\ = 1. 

Here, the quantity p must be regarded as a function of the variable z+ — z_ — 2w, 
where w was defined earlier as a function of p. On the other hand, u is a function 
of — z_ , namely 

= — 7 — — > — 

Observe that the above symmetric formulation does allow the density variable 
to vanish, since the coefficients above remain bounded as the density approaches 
the vacuum. However, since the coefficient in front of dtu in the third equation 
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vanishes with u, we see that the above formulation requires the velocity u to be 
bounded away from the origin which, of course, is not a realistic assumption to put 
on general solutions with vacuum. In Section 4, however, we will discuss a reduction 
of the general initial value problem which ensures this condition after applying a 
well-chosen Lorentz transformation to an arbitrary solution. 

Remark 2.3. The assumed normalization |up = 1 could he relaxed in the formu- 
lation of the system. In fact, if this condition holds at the initial time, then it holds 
for all times, as is clear from the transport equation satisfied by \u\ 

dt\u\ + u- V|u| = 0, 

which follows from (|2.10p multiplied by u. 

3. Properties of Lorentz transformation 

Transformation formulas. We will need to rely on the Lorentz invariance prop- 
erty of the relativistic Euler equations and, therefore, in the present section, we 
collect several technical results about Lorentz transfomations. 

For every U £ K." with U ^ 0, we set U := C//|[/| and we decompose any vector 
X £ M" in a unique way such that 

X = x\\U + Xi_, = a; ■ [/ G R, x^ ■ U ^ Q. 

The Lorentz transformation {t,x) ^ {t',x') associated with the vector U is then 
defined by 

t' ^jiU) (t-e^U-x), 
4 - liU) ix\\ - [/|| t), 
Xj_ — x±, 

where 

7(f/) = 



1 



v/l-e2|J7|2 

is the so-called Lorentz factor. This transformation can be put in an equivalent 
form 

t' = j{U){t-e^U-x), 
3.1 , ^ 

x' = --f{U)Ut +{!+ (7(C/) - 1)U (g) U)x. 

It may be also convenient to use the modified velocity scalar V associated with 
U (following the definition in the previous section) and given by 

T + e|C/|\i/2 



and to rewrite the Lorentz transformation as 

{t' ±eU -x) =e=F^^ {t±eU -x), 
x'^ = x±_ . 

or, equivalently, as 

t' = cosh(eV") t - e sinh(eF) U ■ x, 
(3.2) x' = - sinh(ey) t + e cosh{eV) U ■ x, 

x'^ — x±. 



8 



P.G. LEFLOCH AND S. UKAI 



Recall also that Lorcntz transformations together with spatial rotations form 
the so-called Poincare group of isometries, characterized by the condition that the 
length element of the Minkowski metric is preserved, that is, 

Recall also that the relativistic Euler equations are invariant under Lorentz trans- 
formations which, for instance, can be checked by direct (tedious) calculations or 
from more abstract considerations. The following transformation rule will also be 
useful in the following section. 

Lemma 3.1 (Velocity transformation formula). Let U G M" with U 0, and 
denote by u, u' the fluid velocity vectors in different coordinate systems {t, x), {t' , x') 
related by the Lorentz transformation p.ip . Then, the transformation law for these 
velocity vectors is 

(3 3) - T^hr^ ( - ^ + (^(^)"^ + (1 - ^ « ^) ") 

=: i$(eu,e[/). 

Proof. The fluid velocity vector represents the velocity of a fictitious point-mass 
moving along with the fluid. Thus, u = u{t,x) represents the velocity vector of a 
point-mass located at x at the time t, while u' = u'{t',x') is the velocity of the 
same point-mass in the coordinate system (t' ,x'). Consequently, the vectors u and 
u' are given by 

, ^ , , doc I dx 

Now, in view of (|3.ip , 
, dx' fdt' 



~ ~ d d -1 

= ( - l{U)U + (/ + (7(C/) - l)U ® U) ^) (j{U) {l-,^U- ^)) , 

which, together with (j3.4p . yields (|3.3p . □ 

A technical property on the Lorentz-transformed velocity. In the following, 
we will need to have a lower bound on the fluid velocity vector, so we establish 
here a preliminary estimate. Throughout, e G (0, 1) and all of the constants are 
independent of e. Given tq G (0, 1), we define 

Br, {eu G | e\u\ < tq}. 

Lemma 3.2 (Uniform bounds for the velocity). Given any ro G (0,1) and any 
vector L/ G M'^ satisfying ro < e\U\ < 1, there exist positive constants < 6i < S2 < 
1 depending only on ro and eU , such that the Lorentz transformed velocity p.3p 
has a norm uniformly bounded away from, both, the origin and the light speed, i.e. 

Si < |$(eu,e[/)| < 62 

hold for any eu G Br, ■ 

We observe that the above statement is sharp, in the sense that the constants 
61,62 may approach the endpoints of the interval (0, 1) when tq also approaches 
the endpoints of (0, 1) or when U approaches the endpoints of the interval (rg, 1). 



A SYMMETRIZATION OF THE RELATIVISTIC EULER EQUATIONS 



9 



Proof. To simplify the notation, we introduce the new variables and function 

X:=eu, Z:=eU, WiX, Z) := {eu'l"^ = {^{eu, eU)f . 
It is convenient to choose the coordinate system so that 

Z = (ri,0,0), 0<ro<ri<l. 
Setting X = (Xi, X2, X3) and noting that Z = (Xi, 0, 0), we get 

(3.6) w{x, z) = I YzjTz + - ^(f^)"')^ ®zx-z)^ 

= (T3^ ((^1 - + 7(c/)-^ {XI + xD) . 

The norm \Z\ = ri being fixed, we are going now to compute the extremum values 

of the function W within the domain \X\ < Tq. 
We first compute the minimum by noting that 

Xl - n ^ 2 



An obvious lower bound is W{X, Z) > {n — ro)"^ {1 + rivo) ^. To obtain an optimal 
bound, we compute 



1-nXi (l-riXi)2' 

which is negative for any \Xi\ < ro < ri. Consequently, the minimum is attained 
at Xi = ro, and 

e\u'\=W{X,ZY/^>g{roYI^ 

^^•^^ = =: 5, > 0. 

1 - ron 

Next, to compute the maximum, we set \X\ = r < ro- Noting that 7(C/)~^ = 
1 - = l-rl, we get 

{X,-r^f+'y{U)-\Xl+Xl) 

= (1 - Xiri)2 + x',-l + rl- r\Xl + (1 - r?)(X| + X|) 

= (1 - Xiri)^ + Xl{l - t\) - (1 - t\) + (1 - t\){XI + Xl) 

= (l-Xrrrf-{\-r\){\-\X\^) 

= (l-Xiri)2-(l-r?)(l-r2) 

and, therefore, 

= 1 - /i(Xi,r), /i(Xi,r) = ^ ■ 

Clearly, for each fixed r G (0, ro], the function h{X\^ r) attains its minimum value 
at the point X\ = — r, so that 

(l + rri)2~(l~r2)(l-r2) 



W{Xy,Z) <l-h{-r,r) 



(l + rri)2 
(l + rri)2 - 
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Since 

^ ' (l + rri)3 

for all r > we have 

Finally, by choosing 82 '■— k{roY^^ we obtain the desired inequality and the proof 
of the lemma is completed. □ 

A symmetrization for non-relativistic fluids. For clarity, let us explain our 
strategy to avoid the zelo velocity problem in the simpler case of the non-relativistic 
Euler equations 



dtp + ^d.^^{puk) = 0, 

(3.8) „ 

dt{puj) + ^dx^,{pujUk+pSjk) =0. 

k=l 

Since this is just for the limit case e — > 0, and since the symmetrization 

argument of Section [2] is still valid for this case, setting e = in Proposition 12.21 
yields a symmetrization of p.8p . in the form 

dtz+ + {\u\ + c{p)) u ■ Vz+ + c{p) \u\ iv{E{u) Vu) = 0, 

(3.9) dtz^ + (\u\ - c{p)) u ■ Vz_ - c{p) \u\ tT:{E{u) Vu) = 0, 

2 \u\'^{dtu + u ■ Vu) + c{p) \u\ E{u) Vz+ - c{p) \u\ E{u) Vz_ = 0. 

This symmetrization still has the drawback of having a possibly vanishing coefhcient 
(the velocity) in the third equation. 

Now, recall that the non-relativistic Euler equations (j3.8p are invariant under 
Galilean transformations and introduce the coordinate frame translated at some 
given velocity U . We denote the new variables and unknowns with the symbol jj, 
that is, 

(3.10) —X- Ut, p« := p, J :^ u - U. 
We first obtain, by virtue of the Galilean invariance. 



9V + ^C(p«4)=0, 

(3.11) 

dl{p^u^)+j2dl{p^u]ul+pH,,)^Q. 

k=l 

and then we note that the following symmetrization of (|3.1ip can be deduced from 

the same argument as above: 

(3.12) 

dlz\ + (I J I + c(p")) ■ V»z| + c(p«) \u\ tT{E{u^) V"2») = 0, 

dlzl + (|u«| -c(p«))m- Vzi - c(p«) |w«|tr(i;(M«) V*M«) = 0, 

2 \u^\^{dtu^ + w« • V'u") + c(p«) \u^\E{u^) V»zi + c(p«) \u^\E{u^) V^zi = 0. 
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Now, the advantage of the symmetrization (|3.12|) . in comparison with (13. 9p . is 
obvious: In view of (j3.10|) . never vanishes as long as u remains bounded, provided 
the reference velocity U can be chosen so that, say, 2\u\ < \U\. In turn, Kato's 
theory ensures the local well-posedness for the system p.l2p and, as a consequence, 
for the original system (|3.9p . 

In the next section, we will see that the same strategy works for the relativistic 
case provided Galilean transformations are replaced by Lorentz transformations. 



4. LOCAL-WELL-POSEDNESS THEORY 

Relying on the symmetric form discovered in Proposition 12.21 we are now ready 
to establish the main results of the present paper. We denote here by 
the uniformly local Sobolev space of order r > 0. (Recall that, by definition, the 
Sobolev norm in these spaces is computed on unit balls with arbitrary center varying 
in M".) 

Theorem 4.1 (Local-in-time solutions in Sobolev spaces). Consider the relativis- 
tic Euler equation for an equation of state p — p{p) satisfying the hyperbolicity 
condition (jl.3p together with the following condition near the vacuum 

(4.1) limsup-^^ < oo. 

P^o W{p) 

p>0 

For every constant M > 0, there exists k g (0, 1) such that the following property 
holds. Given at t = an initial data p, u belonging to the Sobolev space iJ^; with 
r > I + n/2 and satisfying the constraints 

< p < Af, lup < K, 

there exists a unique local solution p, u to the corresponding initial-value problem, 
which is defined up to a some maximal time T > and satisfies 

p, u G C([0, T), HliiW)) n ([0, T), Hll\W)) 

and 

p>0, e^\u\^ <l. 
The rest of this section is devoted to the proof of Theorem 14.11 

Step 1. By assumption, the initial velocity scalar is bounded away from the light 
speed and, in consequence, thanks to Lemma [3?^ we can find a vector U £ M" with 
sufficiently large norm \U\ such that the transformed fluid velocity It', defined as in 
(USD by 

u' = l-—(-fiU)-\+{l-j{U)-'){U®U)u-u), 

(4.2) l-e^U-u\ J 

= -<fieu,eU), 
e 

is bounded and bounded away from the origin. Precisely, for some constants < 
Si < 62 < 1 we have 

(4.3) ^1 <e|?l'| <^2. 
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Step 2. Recall that the general theory established by Kato [3j covers symmetric 
hyperbolic systems of the form 

n 

(4.4) AQ{W)dtW ^^Aj{W)djW = Q, 

i=i 

where the (d x d)-matrix fields A[j,Aj are real- valued and symmetric with regular 
coefficients, and the matrix A^ is uniformly positive definite. For the system (|2.8p - 
(|2.10p . we have d = n + 2, W — {z+,z-,u)* (a column vector, the subscript "i" 
standing for transposition), and 

(4.5) 

(oq \ / aiUj a2|M|ej 

6o , A^(Ty) = b2Uj -a2\u\ej 

co|up// \a2|u|e* -a2|u|e* co|upUj/y 

where 



ao = l + e \u\c{p), bo = 1 ~ e \u\c{p), cq = _^ _ ^2|^|2 ' 

oi = j-r^ (|u| + c(p)), bi^ -——-— {\u\^c{p)), 

1 — e''c(p)\u\ 1 + e-^c(/9)|u| 

a2 = c{p), Cj = {E.ji{u),Ej2{u), . . . Ejniu)), 

and we recall that / = {6ij) denotes the n-dimensional identity matrix. Recall that 
we are interested in the initial-value problem associated with (|4.4p where initial 
data are prescribed on the initial hyperplane 

Hq: t = 0. 

First, observe that 

(4.7) {Ao{W)t0^ao\iif + bo\a\' + co\uf\^f, 
where (•, •) denotes the Euclidian inner product in K."+^ and 

e = (ei,6,---,en+2) = (ei,6,0 gR"+^ e = (6,---,e«+2) eM". 

As was already noted before Remark l2.3l and is readily seen directly from (|4.6p . the 
matrix Aq can be positive definite only if the velocity u never vanishes. In other 
words, provided the initial velocity u is bounded away from 0, according to Kato's 
theory, a local-in-time solution exists and is unique in the uniformly local Sobolev 
space for s > f + n/2. As stated in Section 1, however, this lower bound on 
the velocity is not physically realistic. 

On the other hand, from the physical view point, the zero velocity is not a special 
value and, in any case, one should be able to recover the strict positivity property 
for the matrix Ao. At this juncture, recalling the strategy presented at the end of 
the preceding section for the non-relativistic Euler equations, we propose ourselves 
to apply a Lorentz transformation. 

Using the Lorentz invariance property f Lemma 13. ip of the Euler equations, we 
see that the symmetric formulation (|2.8p - (|2.10p can be also expressed in the trans- 
formed coordinates {t',x') defined by p.ip . that is, 

n 

(4.8) Ao{W')dt'W' + J2AAW')d,'W' ^0, 
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where W = (z'^, z'_,u') is defined from the transformed unknowns {p',u'). (Of 
course, the mass density remains unchanged but is now regarded as a function of 

After this transformation, the expression (|4.7p becomes 

(4.9) (AoiW) = «o iCiP + ^'o I6P + c',\uf ICP, 

with aQ,6Q,Co defined by l|4.6p with (p, m) replaced by {p' ,u'). In view of the lower 
and upper bounds (|4.3p . we conclude that the transformed matrix Aa{W') is pos- 
itive definite in the coordinate system {t',x'). Hence, Kato's theory applies to the 
initial value problem for (|4.8p . without any assumption on the fluid velocity, but 
provided initial data are imposed on the initial hypersurface t' = 0. 

In contrast to the non-relativisitic case, however, this is not the end of our 
discussion, since = is not the hypersurface of interest. This is due to the 
fact that, in the relativistic setting, the initial plane Tio is not preserved by the 
transformation (|3.1[) . 

Step 3. In fact, the initial hyperplane is mapped, in the new coordinate system 
(t',a;'), to the "oblique" hyperplane 

Ho : t' = -e^U-x'. 

In order to prove local well-posedness for the oblique initial-value problem to (|4.8|) 
with data prescribed on Tig, it is convenient to introduce a further change of coor- 
dinates 

(4.10) t" ^t' + (:^U -x', x"=x', 
which maps the hyperplane Tig to the hyperplane 

Hg : t" = 0. 
This transformation puts the system (j4.8p into the form 

n 

(4.11) Ba{W") d^>>W" + Bj(W") d^nW" = 0, 
where W"{t" ,x") = W'{t',x') and the new matrix-coefficients are 

n 

Bo{W") - Ao{W') + e'Y^ U, A,{W'), 

B,{W")=A,{W'), j^l,...,n. 

(Note in passing that using a Lorentz transformation instead of (|4.10p would be 
physically more natural, but would lead to precisely the same matrix Bq and slightly 
more complicated expression. 

We are now going to establish that the matrix Bo{W") is positive definite for 
data that have bounded mass density and whose velocity scalar is bounded away 
from the light speed. Provided this is checked, Kato's theory then applies to the 
initial value problem for (|4.1ip on the hyperplane Hq, which is the one of interest. 

At this juncture, it is worth recalling the standard fact that the definite-positivity 
property above implies that the oblique hyperplane Tig is a non-characteristic hy- 
persurface for the hyperbolic system (|4.8p which also is sufficient to imply local 



(4.12) 
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well-posedness. Namely, the matrix Bo{W") can be written as 

n 

Bo{W") =Y,''aAa{W"), 

a=0 

where the vector ly — {vq, . . . , i/n) — (1, e^U) G R"+^ is normal to Hq. Thus, the 
positivity property of Bo{W") implies also that det(i?o(W^")) 7^ 0. 

Let us summarize the expressions we need here. An easy computation with 
(|4?T2| shows that 

{Bo{W") =K + a\e\U ■ Sf'))^? + 2a', ^ |u'|(C/ • 6 
(4.13) + % + fele^ (f^ . ^2 _ 2„^,2 |^,| . (Ei^y^l)) ^2 

where we recall that ^ = (Cii ^2, ■ • • , Cn+2) = (Ci, ^2, G K"^^. The primed quan- 
tities Cq etc. are still defined by (|4.6p but with m replaced by u' determined by the 
Lorentz transformation associated with the reference velocity \J . 
Hence, we can regard the expression 

as a polynomial in ^ and a nonlinear function in eu, €c{p')^ eU. As before, to simplify 
the notation, we introduce 

X-~eu, Y:=ec{p'), Z := eU, 

to deduce 

g(e; X, Y, z) =(4 + a'Mx, z) ■ + 24 z)| (z • z))o) a 

+ {h'^ + z) • z) e2 - "^ci'^W, Z)\ {U ■ {E{^{X, Z))0) 6 

where (after appropriate rescaling in e for the coefficients a^, 6^, a,) 
(4.14) 

a;, = i + l$(x,z)|y, 6;, = i-|<i>(x,z)iF, c[, = ^-^^^, 

(|<&(x,z)| + r), 6; = _^_ZL^(|ci,(x,z)|-r), 



' i-r|$(x,z)p' ^ ' ^' ^' ^ i + y|$(x,z) 

with $(X, Z) := $(X, Z)/|$(X, Z)|. Replacing the coefficients a'Q,a[, etc by their 
values, we finally obtain 

(4.15) Q{C,X,Y,Z) = Q,ei + 02^2' + Q3 lei' + Qi3(6,0 + Q23(6,e), 
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where 



Qi := 1 + mx, Z)\Y + Z) ■ Z) Z)\ + Y), 



Y\^{X,Z)\ 

1 - Z)\Y + Z) . (|<i>(X, Z)l - y), 

2|$(X, Z)p 



Qi3(6, := 2r |<i>(X, Z)\(Z ■ {Ei^{X, Z))t)) ix 

Q23(6,0 := -2r |$(x,z)| (z- 6, 

Recall that e(^(X,Z)) = I -'^{X,Z)® Z). We are interested in the range 
where X, Z G R" have norm bounded away from 1, and Y remains in a bounded 
closed subset of [0, 1). 

It remains to check the following purely algebraic result. 

Lemma 4.2 (Uniform positivity property). For any given Yq G (0, 1), there exist 
r* G (0, 1) and Z e M",r, < \Z\ < I such that 

(4.16) Q{^,^;X,Y,Z)>co\e, 
holds for all Y e[0,Yo], \X\ <r^. 

Proof. In the below, we fix Ya G (0, 1) and Y G [0, Ya] is an arbitrary number. It is 
convenient to choose the coordinate system so that 



ro = \X\, Z = (ri,0,0), < ro < ri < 1. 



Set 



Then 

and Lemma [3.21 savs that 



1 - rori 1 + rori 

1 



< \^{X,Z)\ < $*. 

A simple computation shows that 

z • Hx, z) = Yx^f \^\ = T^'^ =^ ^ < 0' 

1 — A • Z i — Xiri 

which, together with W{X, Z) in (|3.6p . leads to 

R 



S : = Z-<^iX,Z) = 
= \Z\ 



iy(X,Z)l/2 



■((Xi-ri)2+7([/)-2(X| + X|)2)^/' 
Thus, we obtain 

-ri--|Z|<5<0, 
and the equality in the first inequality is realized for X = (^i, 0,0). 
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As a consequence, we have, 

1 - Z)\'Y^ + S{1 - Y^) Z)\+Y) 



Qi = 



l-Y\^X,Z)\ 
(1 - + (1 _ y2) (1 ^ S{\^{X, Z)\ + Y)) 



\-Y\<^{X,Z)\ 

(4.17) > (1 - Y^) (1 - ri($* + Fo)) =: q^. 

Therefore, we see that if 

(4.18) A;o :=ri($*+yo) < 1, 

then gi > and so, Qi > 0. 
Similarly, 

1 - Z)|2f2 + _ y2) (|$(j^^ Z)| - y) 



02 



> 



1 + y|$(x, z)\ 
(1 - |^(x, z)ny2 + (1 - Y^) (1 + 5(|^(x, z)\ - y)) 
i + y|$(x, z)| 

(1 - If) (1 _ niax(0, |<i>(X, Z)| - F)) 

1 + Yo** ' 
Therefore, we have 

(l-yo2)(l-ri^*) _ 

TTi^^J > °' 

since \^[X,Z)\-Y < \^{X,Z)\ < $*. 
On the other hand, clearly we have 

and 

IQ13 (a , f) I < 2y I z) I |z| la 1 1^ < 2yo$Vi 16 1 1^ , 
i<323(a,oi < 2y i$(x,z)iiz|iaiia < 2io$Viiaiia- 

Set now 

q4 := ^0**7-1 

and define a quadratic formula of three variables (a;, y, z) G M'^, 
Q*(x, y, z) = qix^ + q2y^ + q^z^ - 2q4{x + y)z. 

Then, we obtain 

Q(ax,F,z)>g4iai,i6i,ici) 

for any ^ G M"+2; that is, Q is positive definite if so is Q*. 
Since gi, (72 > 0, we can write, for any k € (0, 1), 

Q*{x,y, z) = K{qix'^ + q2y'^) + (1 - K)qi (x - z) 

V (1 - K)qi / 

+ (1 - My - TT^^y + - - 

\ {1-K)q2 / V (1-K;)g'i (1-K)g2>' 



z', 
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and since k G (0, 1) is arbitrary, we can conclude that Q„ is positive definite if and 
only if 

(4.19) D,:^qs-^^^>0. 

Let G (0,1) be a number to be determined later. Let a £ (1,1/r,) and set 
ri = ar, . Then, observe that for any rp G [0,r,], 

n-rp ^ (a - l)r^ + ri (a + l)r^ ^ ^ 

= -j > -j —, $ = -r- <— < a+l r*. 

1 — rpri 1 - arj 1 + rp^i 1 + rpri 



Hence, 



and 



ko = ri($* + Yo) < ar^{{a+l)r^+Yo) < a(a + 2)r, fcir,, 
gi = (1 - Y^)il - ko) > (1 - ro')(l - a(a + 2)r,) =: ifi(r*). 



1 + Yon 2 



1 - $2 V / (1 - ar2)2 - (a - l)2r2 

94 = ^o***^! < i"o(a + l)arl =: K^rl 

Note that ki , K4 > are independent of and that all the above inequalities hold 
for all To G [0,r*]. Observe that if r, is sufficiently small, then A'i(r*), i4r2(r,), 
K3{rf) are positive and 



(4.20) > (^Ksin) 



holds for all rp. It is easy to see that when — > 0, we have 

fcp^o, ifi(rO i-rp2, 

and 

if2(r0^^(l->^'), K^{r.)^2{a-lf. 

We can now conclude from (|4.20|) that for aJi a > 1, there exists G (0, 1) such 
that or* < 1, fcp < 1, I?* > 0, which completes the proof of Lemma 14.21 with 
ri = \Z\ = ar<,. □ 

In turn, Kato's theory guarantees the existence of a solution defined in a small 
neighborhood of this hyperplane Tig . Making the transformation back to the origi- 
nal variables, we obtain a solution in a small neighborhood of the initial line t = Q. 
This completes the proof of Theorem 14.11 



Step 4- We need to show that the density p remains non-negative. In fact, since 
the solution remains smooth enough, say of class in space and continuous in 
time, we can define the characteristic curves by the Cauchy-Lipschitz theorem 



y{t)^u{t,y{t)), t>0. 
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Hence, by writing the w-equation as 
[dtw + u ■ Vw) 



and integrating along the characteristics we obtai 




^(wit,yit))^0{l)c{pit,yit)). 



Here, the (bounded) quantity 0(1) involves only the sup-norm of first-order deriva- 
tives of the solution. 

Now, in view of our assumption (j4.ip we have 



on any compact set p G [0, R] and for some constant Cr. So, we deduce that 



along every characteristic, which, in particular, shows that w remains non-negative. 
This completes the proof of Theorem 14.11 



Following [7] we define a concept of solutions, in which the velocity vector is 
required to satisfy a (non-degenerate) evolution equation even in the presence of 
vacuum, as follows. 

Definition 5.1 (Notion of tame solution). A measurable map {p,u) : [Q,T] x 
[0, cxd) X [— e^^,e^] is called a tame solution of the relativistic Euler equations if 

• (p, u) is a solution of class of the Euler equations, 

• w is also of class , and 

• the equation dtu + u ■ Vu = holds in the interior of the set {p = O}. 

Relying on this definition, wc can establish that the support of a solution does 
not expand in time. 

Theorem 5.2 (Property of the support of a tame solution). Consider the rela- 
tivistic Euler equation for an equation of state p — p{p) satisfying the hyperbolicity 
condition (|1.3p together with the following vacuum condition (14. ip . If {p,u) is a 
tame solution of the relativistic Euler equations and has compact support, then its 
support does not expand in time, that is, 

supp{p,u){t) C supp{p,u), t £ [0,T]. 

Proof. We now consider the relativistic Euler equations in the form 



c{p) < Crw{p) 



w{t,y{t)) 



u;(0,y(0))e' 



,to(i) 



5. Support of tame solutions 



n 



BoiV) dtV + BjiV) djV = 



i=i 

where the map V := {w,v,u) is of class and satisfies 

\dtV\ <Ci snp\BoiV)-^B,iV)\. 



j 
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Here, the constant Ci depends on the sup norm of first-derivatives of the solution. 
In view of the explicit expressions of the matrices Bq and Bj (see Section [2]) we 
obtain 

snp\Bo{V)-'Bj{V)\<C2\V\, 
j 

where we have used our assumption (|4.ip . Now, by Gronwall's lemma we get 

\V{t,x)\ <e'^*\V{0,x)\, t>0,xeW\ 
which completes the proof of the theorem. □ 
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